High temperature anomaly of the conductance of a tunnel junction 
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The linear conductance of a tunnel junction in series with 
an ohmic resistor is determined in the high temperature limit. 
The tunneling current is treated nonperturbatively by means 
of path integral techniques. Due to quantum effects the con- 
ductance is smaller than the classical series conductance. The 
reduction factor is found to be nonanalytic in the environmen- 
tal resistance for vanishing resistance. This behavior is a high 
temperature manifestation of the Coulomb blockade effect. 

PACS numbers: 73.23.Hk, 73.40. Gk, 73.40.Rw 



In the last decade or so tremendous progress was made 
in our understanding and control of Coulomb blockade 
phenomena |3,0|. These effects arise in nanostructures 
with small capacitances C when the single electron charg- 
ing energy Ec = e^/2C becomes significant. So far work 
has concentrated on the low-temperature, typically mil- 
likelvin region to satisfy the seemingly obvious condi- 
tion ksT <g; Ec for the manifestation of charging ef- 
fects. However, recently Pekola and coworkers [pi have 
demonstrated that there are remarkable signatures of the 
Coulomb blockade effect also for temperatures in the few 
Kelvin region well above the charging energy. In fact, 
charging phenomena in the high-temperature region can 
be used for precision thermometry g. 

Further, in most of the theoretical studies of charging 
effects in metallic nanostructures it is assumed that the 
tunneling resistance Rt of the tunnel junctions in the 
device exceeds the von-Klitzing resistance Rk = h/e^ « 
25.8fcr2. Under the condition Rt » Rk quantum fluctu- 
ations of the charge Q are strongly suppressed. Hence, in 
a sense this latter condition ensures that charging effects 
are not washed out by quantum fluctuations in much the 
same way as the above-mentioned temperature constraint 
ensures stability against thermal fluctuations. 

It is now interesting to examine whether signatures of 
the Coulomb blockade effect persist when thermal ffuctu- 
ations are large and the tunneling resistance may become 
small. This question will be addressed for the fundamen- 
tal problem of an ultrasmall tunnel junction biased by 
a voltage source. The single junction system was stud- 
ied previously ^ in the limit Rt ^ Rk where tunnel- 
ing can be treated perturbatively. While the theoreti- 
cal approach employed in these earlier studies does not 
readily extend to the case of strong tunneling, we know 
from these investigations that the impedance of the leads 
connecting the tunnel junction with the voltage source 
crucially affects the measurable conductance. In fact. 



charging effects in single junctions only arise from the 
coupling of the tunneling electrons to the modes of the 
electromagnetic environment of the junction ||6|. In the 
model studied in this work the environmental impedance 
will be assumed to be ohmic. 
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FIG. 1. Circuit diagram of a tunnel junction with ca- 
pacitance C and tunneling resistance Rt biased by a voltage 
source V via leads of impedance R. 

Specifically, we consider the circuit depicted in Fig. [|. 
We aim at the zero-bias differential conductance G = 
dI/dV\v=o where / is the current and V the appHed 
voltage. Our approach is based on the path integral for- 
mulation of the tunnel junction H and employs the gen- 
erating functional 



Z[£^] ~ tr T-r exp ■ 



up 



dr [H - /e(r)] 



(1) 



where H is the Hamiltonian of the system for V = and 
/ the microscopic current operator for charge flow in the 
leads. Tr is the time ordering operator for the imaginary- 
time variable r. Introducing a phase variable ip m which 
is related to the voltage across the tunnel junction via 
a Josephson-type relation, the generating functional (U^) 
may be written as a path integral 
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D [if] exp 
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S[V,^] 



with the effective Euchdean action 

s[ip,C\ = ScM + STM + SR[^,i] 

where 



(2) 
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describes Coulomb charging of the junction and 


St['A = 2 / dT dr' a{T - t') sin^ 


Xt) ~ip{T')' 
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(4) 

quasi-particle tunneUng across the junction M]. The ker- 
nel a{T) is determined by the tunneling resistance Rt 
and may be written as 
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where the Un — 2-Kn/fij3 are Matsubara frequencies and 
The last term in the action (pf) reads 

Sr[V.S] = -J dr j dT'k{T~T') 



v{r) + -^ar)-vir')-^aT') 



(5) 



For ^(r) =0 this is the well-known influence functional 
of an Ohmic lead resistance R [BJ with the kernel 
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where 



k{vn) 



h Rk 



Note that the influence functional (|^) of the lead resis- 
tance is at most quadratic in the phase while the tun- 
neling term (0) with a trigonometric phase dependence 
renders the path integral (||) non-Gaussian. This is due 
to the discrete charge transfer across the tunnel junction. 
The conductance G can be calculated from the Kubo 
formula 



G = lim - — Im 



lim , / dTe*'^""(/(T)/(0)} 



(7) 



where the current-current correlator (/(t)/(0)) reads in 
terms of the generating functional (M 



(/(r)/(0)) 



S'Z[^] 



z[o] saT)sm 



?(t)=0 



Now, using the path integral representation 
Eq. (P), the correlator may be written as 



(/(r)/(0)) = | /i^Mcxpl-i^M 



(8) 
and 

(9) 



X {2jk{T)+I[^,T]I[ip,0] 



where S[(p] — S[ip,£, =0] and 



(10) 



is the partition function. Further, the current functional 
I[ip, t] is given by 



2e /■^'^ , 
I['P,t] = — dr' k{T-T')(p{T'). 



(11) 



To calculate the conductance from Eq. (|^), we first need 
to determine the Fourier components 



G{v^)^ / dTe'--^I{T)I{Q)) 
Jo 



(12) 



Using Eqs. (|), (g), and (y), one finds that C(i^„) 
Gi{iyn) + C2{vn) where 



Ci(i/„)=2-fcK) 



(13) 



and 



C2{yn) 



|yi?Mexp|-i5M|F[^,z.„]. (14) 



Here the functional 

4e2/3 
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F[ip,l'n] = ^—k{Vn)ip{Vn) ^ k{Vrn)'f{Vm) (15) 



m— — oo 



is given in terms of the Fourier coefflcients ^(f„) of the 
phase variable. Since the action S[ip\ is invariant under 
a global phase shift, we may put ^(0) = 0. 



According to the decomposition of C{vn) into (|13[) and 
(|lj) , the conductance (Q) splits into G = Gi + G2 where 

Gi = lim -^ Im %-k{-iuj + 5) = -. (16) 

i^^o fiLo h R 

Since the kernel fc(T) has the analytic properties of an 
equilibrium correlation function, the analytic continua- 
tion of k(yn) is unique H. 



To proceed we first note that in terms of the Fourier 
coefficients (p{i'n) the action functional may be written 



S[ip] = So[(p] + ^S2k['f] 

k=2 

where 

oo 

So[ip\ = n^\{Vn)\'^{l'n)\'^ 
n=l 

with the eigenvalues 



(17) 
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Here g — Rk {Rt^ + R^^) is the dimensionless parallel 
conductance of the tunneling and lead resistances. Fur- 
ther, 



2/c-l 



^'■[''i^w^C"^'-"'''' 



1=1 



J2 a(~I]p^l'^np)'^(^ni)--' 
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Now, in the high-temperature limit we may expand about 
the Gaussian action So[lp] using 

eM~S[^]/h) = eM~SQ[^]/h)[l - S4^]/h (19) 

-Se[^]/n - Ss['f]/n + S4[f]V2h^ + •••]. 

The ratio of path integrals (|4|) can then be evaluated by 
means of 

J D[:f]exp{-j-SQ[ip]} 'f{i^k)(pM _ Sk-i 



/i5Mexp{-i5o[^]} 



Xii^k) 



and Wick's theorem for higher order products of the 
Fourier coefficients tpiyn)- Since 
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PEc 



\{vk) 27r^k^ + gl3Ec\k\ 

is of order PEc, the expansion (Km gives a high- 
temperature series in powers of (3Ec ■ 

Proceeding along these lines, one obtains from 
Eqs. (|l|) and (|l|) 






a{iy„+„i) - a{vn) - a{vm) 
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In view of Eqs. (Q) and (p^), the high-temperature se- 
ries for C2{vn) yields a corresponding expansion of the 
conductance 



G2 



lim - — ImC2(— itj + 5) . 



(21) 



However, the analytic continuation and the limit w — > 
involve a subtlety. As is readily seen from Eq. (p^), a 
factor l/\{vn) in Eq. ( po| ) gives in this limit a factor 
2TTi/gPfiuj. While the l/w-divergences are canceled by 
tj-factors stemming from the numerators in Eq. (pOj), as 
a net result, each l/A(i/„) factor reduces the order in 
(3Ec of the corresponding term in the high-temperature 
expansion of G2 by one order. Hence, the quantum cor- 
rection to G2 of a given order in PEc depends on higher- 
order terms in the series expansion of C2{vn)- 

An analysis of contributions of all orders |lO| shows 
that a term of the series (n9[) with a product 

S2k.Af]S2k2[f] ■ ■ ■ S2kA'P] 

of quartic or higher-order actions gives quantum correc- 
tions to G2 of order (/3£;(^)'=i+'=2+- +fe<!-^ g^^^ qJ higher 

orders. This proves that the terms of the expansion of 
C2{i^n) given explicitly in Eq. ( po|) suffice to calculate the 
leading order quantum corrections to G2. 

Now, combining Eqs. (|2y) and (^Ij) one finds after some 
algebra the explicit high-temperature result for G2 that 
may be added to Eq. (16). This gives for the conductance 



G = 



1 



Rt + R 
Rt 



1 



7 -I- ^"(1 + u) 



Rt + R 

' 0{l3Ecf 



+ V/(i 



pEc 



where 7 is Euler's constant and 

_g(3Ec _ Rk{Rt + R) f3Ec 



27r2 



RtR 



27r2 



(22) 
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Further, ^p{z) and ip'{z) denote the digamma function 
and its derivative, respectively. 

The predicted quantum corrections to the classical con- 
ductance l/(i?T + R) reduce the conductance due to the 
Coulomb blockade effect. In Fig. g the result (|2^) is de- 
picted for fixed temperature and tunneling resistance Rt 
as a function of R/ Rk- Typical lead resistances are in 
the range of lOOfi so that R/ Rk is usually very small. 
Then, the parameter u defined in (p3) is large even for 
temperatures T well above Ec/ks, and one can use the 
asymptotic expansion of the digamma function to give 
for R<€.Rt,Rk 



G 
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1 + 2 j^ In 
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\Rk 
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hi{PEc)\R 
Rk 



(24) 



This shows that for small lead resistance the influence of 
the electromagnetic environment cannot be treated as a 
weak perturbation since a finite resistance gives rise to a 
reduction of G that is non-analytic in R. This large effect 
of even a small lead resistance is clearly seen in Fig. 0. 
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even in the high-temperature range where ksT is large 
compared to the single electron charging energy Ec- In 
particular, in the experimentally relevant range of small 
environmental impedance R, the corrections to Ohm's 
law are significant, since G{R) is nonanalytic at i? = 0. 
In view of the fact that the conductance can be measured 
very accurately, the predicted effects should be readily 
observable. 

The nonclassical behavior of the conductance is also 
important for the precise determination of junction pa- 
rameters from data in the "classical" regime. In this con- 
text it is important to note that the result ( p2| ) remains 
valid for small tunneling resistance Rt ^ Rk as long as 
NRt 3> Rk where N is the number of transport chan- 
nels. Since N is typically above 10^ for metallic tunnel 
junctions, the theory can be used to determine junction 
parameters also in the region of strong tunneling. 
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FIG. 2. The ratio of the conductance G and the classical 
conductance 1/{Rt + R) is shown vs the environmental resis- 
tance R in units of Rk for Rt = Rk and pEc = 1/4. The 
full line shows the result (221) , and the dotted line depicts the 
approximation (M) for moderate-to-high resistance. 

On the other hand, for values of R/ Rk of order 1 or 
larger, u becomes small in the high-temperature limit 
and the digamma functions in Eq. (p3) may be expanded 
to yield 



G 



1 



Rt + R 



Ri 



PEc 



Rt + R 3 



0[il3Ecf,u(3Ec 



We note that in the limit of weak tunneling, Rt ^ 
R,Rk, this result may be approximated further to read 
G = i?^^[l - (3Ec/'i + 0{f3Ecf], which is in accordance 
with the conventional approach P,o| based on perturba- 
tion theory in the tunneling term. Within this framework 
higher order terms in PEc can be evaluated explicitly 
pT[ . However, the conventional approach misses the fact 
that the quantum corrections to Ohm's law vanish in the 
limit R/Rt -^ oo according to (p5[). For large R many 
charges fluctuate across the tunnel junction on the rele- 
vant time scale RG and lead to a suppression of charging 
effects. 

In summary, we have shown that the conductance G of 
a tunnel junction is affected by Coulomb charging effects 
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